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ABSTRACT 

Flexion is a non-linear gravitational lensing effect that arises from gradients in the con- 
vergence and shear across an image. We derive a formalism that describes non-Hnear 
gravitational lensing by a circularly symmetric lens in the thin-lens approximation. 
This provides us with relatively simple expressions for first- and second-fiexion in 
terms of only the surface density and projected mass distribution of the lens. We give 
details of exact lens models, in particular providing fiexion calculations for a Sersic-law 
profile, which has become increasingly popular over recent years. We further provide 
a single resource for the analytic forms of convergence, shear, first- and second-fiexion 
for the following mass distributions: a point mass, singular isothermal sphere (SIS); 
Navarro-Frenk- White (NEW) profile; Sersic-law profile. We quantitatively compare 
these mass distributions and show that the convergence and first-fiexion are better 
indicators of the Sersic shape parameter, while for the concentration of NFW profiles 
the shear and second-fiexion terms are preferred. 
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1 INTRODUCTION 

Quantifying the exact shape of baryonic and dark mat- 
ter density profiles that form in expanding ACDM cos- 
mologies is an ongoing issue. N-body simulations suggest 
that CDM halos are well fitte d by either Navarro -Frenk- 
White (NFW) density profiles ijNavarro et al.lll997 ) or de - 
pr ojected Sersic-like laws in the form of Einastc lll965ll 
or Prueniel & Simien (1997") density profiles (Navarro et al,l 
1.2004; Merritt et al. 2005, 2006; Graham et al. 2006a,b). Nu- 
merous observations of early-type galaxies suggest their 
luminosity profiles, and hence stellar mass distributions, 
follow either Ser s ic, core-Sersic or Nuker-law models (e.g. 
Lauer et al.''l995'. 2005'; iGraham et"aDl2003l : iFerrarese et al.l 
20061 : Cote et al. 2006,, and references therein), while gravi- 
tational lensing observations have suggested the total mass 
distribution (i.e. baryons plus dark matter) is consis- 
tentl y descr i bed by isothermal s pheres ijTreu fc Koopmanp 
200l 12004 i Rusin et al.' 2003; Rusin fc Kochanekl l2005l: 



been no less intense; only r e cently is the NFW 
file (e.g. IC arlberg et al '! '199?'; ^- 



'van der Marel et al.^ ''20001; 
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Hansen et al.l 


2005 


, '"11 1 T 

: Lokas et al. 


l2006ir iRines & Diaferiol 


l2006l: Woitak et all 


2OO7I: Okabe et al.ll2009l) being favoured 


over the isothermal sphere (e. 


g. Athreya et al.l I2OO2I: 


lEttori et al.l I2OO2I: iKatgert et al.l |2004|). As gravitational 



Kooproans et al., ,2006 
iDve et al.1 \200i 



200i 



.Gavazzi et al., ,20071: Czoske et al.l 



Tu et al.l I2OO9I ). Meanwhile, de 



bate about the mass distributions of galaxy clusters has 
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lensing traces total projected mass, it is an extremely pow- 
erful tool for determining the mass distributions of these 
systems without having to make assumptions about the dy- 
namics o r constitution of the lensing objects (fo r recent re- 
views see ISchneiderir2005l : iHoekstra fc Jainll2008l ). 

Traditionally, the study of weak lensing has been lim- 
ited to linear effects; convergence and shear. These fields 
have the effect that an elliptically shaped source galaxy 
gets mapped to an elliptical image. Therefore, to deter- 
mine information about the lensing object using first-order 
quantities, assumptions must be made about the intrin- 
sic ellipticity of the source galaxy, or a large number of 
source galaxies must be utilised to ensure statistica lly rea- 
sonab l e results can be inferred (see for e.g. Hoekstr a et al.l 
I2OO4I : iMandelbaum et al.l I2OO6I ). Recently however, vari- 
ous authors have begun to consider higher-order lens- 
ing effects known as flexion dGoldberg fc NatarajanI l2002l : 
iGoldberg fc Bacon|[2005l : lBacon et al.ll2006l 1H Flexion comes 



llrwin fc Shmakoval l|2005l . [2006l ') also consider higher-order lens- 
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in two flavours that correspond to various spatial deriva- 
tives of the shear and convergence, implying they are due 
to gradients in the first-order fields across the extent of the 
source/image. Physically, one can think of first-flexion as a 
shift in the centroid of the image with respect to the source 
and se cond-fiexion as crea ting an arc-like structure in the 
image ijBacon et al.ll2006l l. That is, with the inclusion of 
flexion, an elliptical source galaxy gets mapped to a "jelly- 
bean" shaped image for a circularly symmetric lens. In this 
way, flexion provides a better observable than the flrst-order 
flelds as only one reasonable assumption about the source 
galaxy is required - i.e. galaxies are not intrinsically flexe(J!|. 

The flexion of a lensed image is formally calculated 
using multipole moments jGoldberg fc NataraianI I2OO2I : 
iGoldberg fc Leonard! [2OO7I : lOkura et al.l l2007l . I2OO8I I. how- 
ever in this work we choose to treat the gravitational lensing 
variables as fleld variables. That is, we determine the amount 
of convergence, shear and flexion one would measure as a 
function of the distance from the centre of the lensing mass 
and the angle in the sky, ignoring the overall shape and size 
of the source. Determining the change in shape between the 
source/image pair requires a more detailed mapping that 
calculates small changes in the position of numerous light 
rays in the image and source planes. This is a somewhat 
more difficult task that requires numerical methods, which 
is beyond the scope of the present work. However, treating 
the gravitational lensing terms as fleld variables is extremely 
useful and the benefits it purveys are three-fold: 

(i) It allows us to determine differences in the global grav- 
itational lensing properties from various realistic density dis- 
tributions (in particular see figure [2] below). 

(ii) We can determine the relevant lensing terms (conver- 
gence, shear, first- or second-flexion) for discerning between 
individual shapes of density proflles. For example in sec- 
tion S] we show that convergence and first-fiexion are good 
indicators of the Sersic shape parameter, whereas the con- 
centration of NFW profiles can be determined by looking at 
the shear and second-fiexion. 

(iii) The derivation of analytic solutions is a critical first 
step towards studying fiexion through multiple lens planes 
with arbitrary mass distributions. 

This paper is set out as follows; In section[2]we system- 
atically develop the two-dimensional thin-lens gravitational 
lens equation for an arbitrary, circularly symmetric matter 
distribution, deriving the first-order terms in section [T2] and 
the fiexion terms in section \2l3\ In section [3] we consider ex- 
act forms of the matter distributions, writing down analytic 
expressions for the convergence, shear and fiexion for a point 
mass, singular isothermal sphere (SIS), NFW and Sersic pro- 
files in sections 13. 1113. 2113. 31 and 13. 4l resDectivelv. In this way 
we are providing a single resource where the analytic forms 
for the convergence, shear, first- and second-flexion can be 



ing, which they term sextupol e lensing with com ponents sextupole, 
cardioid and displacement. In llrwin et al.l ll2007l ) the authors show 
that sextupole is equivalent to second-flexion and a combination 
of the cardioid and displacement terms is equivalent to first- 
flexion. 

^ Galaxies are not intrinsically flexed providing they are dynam- 
ically relaxed. 



found for a range of useful density proflles. Finally, in sec- 
tion|4]we compare the gravitational lensing effects of each of 
these profiles. We find that the shear and second-fiexion for 
a Sersic-law profile are systematically greater than for NFW 
and SIS profiles, whereas the convergence and first-fiexion 
of each of the profiles are comparable. We further show that 
the convergence and first-fiexion provide excellent tracers for 
the Sersic shape parameter, whereas the shear and second- 
fiexion are better indicators of the concentration parameter 
for the NFW profile. We make some concluding remarks in 
section \5\ 



2 ANALYTIC LENSING FORMALISM 
2.1 The Thin-Lens Equation 

Given that fiexion considers finite source sizes, we require a 
two-dimensional version of the thin-lens gravitational lens 
equation, whereby the mapping between a point on the 
source plane and the image plane are explicitly expressed. 
In this way the thin-lens gravitational lens equation is ex- 
pressed as 

= - -DiSQi, (1) 

where is the impact parameter on the image plane, rji 
is the distance between the origin of the coordinate system 
and the source on the source plane, ai is the defiection angle 
and Ds, Dl and Dls are the angular diameter distances 
from the observer to the source, the observer to the lens and 
the lens to the source respectively (the lens configuration is 
shown in figure Equation ^ can be put into a neater 
form by using angular coordinates, /3i = rji/Ds and 9i = 
Dl, such that 

(i^=e^~ Q„ (2) 

where ai = ctiDLS / Ds is the scaled defiection angle. Equa- 
tion ^ can be re-expressed as a coordinate map between 
the two angular coordinate systems allowing the mapping 
to be expressed as a linear transformation 

P.^A.jOj, (3) 

where Aij :— d(3i/d6j is the Jacobian of transformation and 
summation is assumed over repeated indices. Equation ([3|) 
implicitly assumes there are no gradients in the components 
of Aij across the image. Physically, this implies that the 
convergence and shear are constant across the image, which 
is not necessarily a good approximation when the source 
has a finite spatial extent. This is a reasonable assumption 
if sources are assumed to be point-like, however for realistic 
cosmology one would like to be able to probe deviations of 
the sh ear and convergence across images. IColdberg fc BaconI 
lj2005l ) therefore generalized equation ([3]) to include higher- 
order terms, coined flexion. These fiexion terms account for 
an image's "arciness", and are therefore relevant if one is 
considering sources of finite spatial extent. The non-linear 
expansion of the thin-lens equation is 

Pi ~ AijOj + —DijkOjOk, (4) 
where 
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Figure 1. Lens configuration. Rattier ttian use ttie angles, 8i and ft, to describe ttie position of the image and the source on the lens 
plane respectively, we use the two-dimensional radial coordinates denoted £,i = D^di and rn = Ds/3i. The distances between the source 
plane, lens plane and observer are all measured in terms of angular diameter distances, and a; is the deflection angle. 



D 



ijk 



(5) 



is the gradient of the Jacobian. The components of the Dijk 
matrix make up the components of the first- and second- 
flexion, a point that we discuss in considerably more detail 
in section [231 

We are interested in explicit expressions for the first- 
and second-order lensing terms with respect to positions on 
the image and source planes. It is therefore instructive for us 
to work in the Cartesian (riij^i) system of coordinates rather 
than the angular {f3i,6i) coordinates (see figure [T]). In these 
coordinates, the Jacobian and its gradient are expressed as 



A - - — 



In this article, we restrict our attention to distributions 
of matter which are circularly symmetric when projected on 
to the lens plane. For such a profile, one can express the 
defiection angle as (see for e.g. ISchneiderll2006l l 

4GM(|e|), 



Oi = 



(7) 



where = Vsi+Cf ^^'^ the projected mass, M, is defined 
as the area integral of the surface density, E(|^|), which for 
circular symmetry is given by 



Mm = 2n / E(C')C'dC'. 



(8) 



Equation (UJ can now be expressed in component form, also 
substituting equation yielding 



1 M(iei) 



ttEc 



(9) 



where the critical surface density has been defined according 
to 



Ds 



47rG DlDls 



(10) 



2.2 First-Order Lensing 



The first-order Jacobian, Aij, can now be calculated by dif- 
ferentiating equation ([9]) which, after some algebra is 



A 



where for convenience we have defined the function 

Q(C) :=^E(C)C'-M(C). 



(11) 
(12) 
(13) 

(14) 



The convergence and two components of the shear are 
expressed in terms of the components of the Jacobian as 
follows 



71 = -{All - A22)/2, 

72 = —Ai2, 

K = 1 - (All -f A22)/2. 



(15) 
(16) 
(17) 
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Evaluating the convergence by substituting equations l|lip 
and l(T2ll into ifTfll implies 



m\) 



(18) 



which is the fa miliar version of this equation (see for e.g. 
ISchneiderl200^ 'l. Evaluating the two components of the shear 
using equations I|lllll3p implies 



71 = 



72 



26(1^1) 



ei6 



(19) 
(20) 



To compare the equations for the shear with a more familiar 
form, we convert the above equations from Cartesian coor- 
dinates, (^1,^2), into polar coordinates on the lens plane, 
{R,4>), defined by — Rcostf) and ^2 = RsirKj). This im- 
plies that the convergence and total shear, 7 = 71 + 172, 
are 



7 = 



y ' 
QjR) 

TrEcr-R^ 



g2«'^> _ lilexp {2i(f 



(21) 
(22) 



This is consistent with showing that the shear is a spin-two 
quantity, whereas the convergence is spin-zero. 

Equations l|2ip and l|22p give the convergence and shear 
distributions as a function of the mass profile of the galaxy. 
Whilst the convergence is only dependent on the surface 
density, the shear is also a function of the mass. This im- 
mediately gives the expected result for a point mass lens 
that the convergence is everywhere zero (except at the ori- 
gin), whereas the shear is non-zero and goes as M/R^ (see 
section l3JJ. 



2.3 Flexion 

The real interest in this article lies in the fiexion terms, which 
are given by the second-order Taylor series expansion, equa- 
tion (U}. The three-tensor, Dtjk, is expressed in terms of the 
gradient of the linear Jacobian of transformation in equation 
(IgJ that, given equations l|lH13p . allows us to calculate the 
individual components of the three-tensor in terms of the 
coordinates ^1 and ^2. After much algebra, one finds 



D2 



D2 



D2 



2Dl^i (e 



ttE^^IC!'' 



Q 



7rE,,|C|6 
-2Dl^i (Cf - 



Q 



TrEcrlCI'^ 

-2Dl^2 (3C? 



7rE,,lC|6 



s 



2Dl^I dE 
2DlC?6 rfE 
2Dl^i^^ rfE 

Scrimp d\^\ 

2DL^i rfE 



(23) 
(24) 
(25) 
(26) 



where for the remainder of the article E — E(|^|), Q = 
Q(|5|) and M = M(\(,\) unless otherwise explicitly stated. 
As shown bv lBacon et al.l lj2006l ). the three-tensor, Dijk, can 
be expressed as the sum of two other tensors, Dijk ~ J^ijk + 
Qijk, which can be written component-wise as 



3J"i 
T2 



T2 



(27) 



Tij2 

Giji 

Gij2 



T2 

Gi 
Q2 

Q2 

-Qi 



Tx 

3.^2 

Q2 

-Gi 

-Gi 
-G2 



(28) 
(29) 
(30) 



In the above, J- = J-i + iJ-2 and G = Gi + iG2 are know n 
as first- and second-fiexion respectively ijSacon et al.ll20o3 ). 
First-flexion is a spin-one quantity that measures the shift in 
the centroid of the image, and second-flexion is a spin-three 
quantity measuring the "arciness" of the image. 

Inverting the above system of equations implies we can 
express the components of the flrst- and second-flexion in 
terms of the components of Dijk', 



(31) 
(32) 
(33) 
(34) 



:Fi = 


-i(I>lll +I>22l), 


J^2 = 


-i(I>211 +I>222), 


Gi = 


-^(^111 -3D221] 


G2 = 


-i(3D2ii -1^222) 



Therefore, utilizing equations I|23II26II . II31II and 
show that the components of the first-flexion are 

Dl 

Ecr 9^1 ' 

Dl aE 

Ecr 9^2 



:F2 = 



(35) 
(36) 



That is, the first and second components of first-flexion (i.e., 
J^i and J^2 respectively) are the directional derivatives of 
the surface density. Note that the flrst-flexion terms do not 
include any functions of the mass, which is again consistent 
with the flexion due to a point mass being zero (see section 
13. ip . The total flrst-flexion, T = T\-\- iT2, is given by 

Dl dE 



Dl dT,{R) 



(6+^6) 



i0 



Ecr dR ^ ^ 

imply ing the flrst-flexion is a spin-one fleld ijBacon et al.l 
l2006( l which is the gradient of the surface density. Given 
that the convergence is proportional to the surface density, 
this is an equivalent way of saying that the flrst-flexion is 
the gradient of the convergence. 

Following the same procedure outlined above, one can 
show that the components of the second-flexion are 



Gi = 



G2 = 



TTEcrl^i^ 

Dl^2 (3C? ~ gg) 

TTEcri?!'^ 



dE 3 
d^ iti3 



4Q 



4Q 



(38) 
(39) 



These form the total second-flexion term, G ~ Gi + iG2, as 
^ Dl ' 



TrEcrjCr 



rfE 3 



4Q 



(?i + *6)' 



(40) 



Finally, expressing the second-flexion in polar coordinates 
and expanding Q to show the explicit dependence on the 
surface mass density and the projected mass one flnds 
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dE{R) 4E(i?) ^ AM{R) 



dR 



R 



TVR3 



(41) 



This expression is relatively simple and shows the depen- 
dence of second-flexion on the surface density, its diver- 
gence and the projected mass distribution. Moreover, one 
can explicitl y see t hat this is a spin-three field, as shown in 
iBacon et all l|2006h . 

It is pertinent to note that ISchneider fc"Erl lj2008l l 
showed that first- and second-fiexion are not observable 
quantities due to the mass-sheet degeneracy. Instead, they 
showed that the observable spin-one and spin-three fields 
are the reduced flexion terms; 



1 - K 



and 



(42) 



where g = "//(I — k) is the reduced shear and JT* = J^i— iT^ 
is the complex conjugate of first-fiexion. Whilst it is ben- 
eficial to keep the mass-sheet degeneracy in mind, a ma- 
jority of work on mea suring fiexion i n real images us- 
ing e i ther the shapelets (Refregierl 12003 : iRefres ?ier fc BaconI 



2002 



2007 



Goldberg fc Bacon 2005.1 or HO LICs ijOkura et al.. 
20081 : iGoldberg fc Leonardl 120071 1 approaches have 



been based on decomposing the higher-order components 
of the images into T and Q. Therefore, for this article we 
continue to focus on the fiexion quantities, T and C/, rather 
than the terms expressed in II42II . 



2.4 A Note On Units 

The system of coordinates we are using in this article are 
somewhat unconventional in the sense that they use dis- 
tance, rather than angular coordinates on the image and 
source planes. We do this as we believe these units are more 
conducive to numerical applications, and are also more de- 
scriptive to the reader. However, our approach implies it is 
worth spending some time elucidating the dimensions that 
we are using for each of the derived quantities. 

The surface density is defined as the projection of the 
three-dimensional density distribution, implying it has units 
of mass per unit area. Integrating this according to equation 
(HI implies the two-dimensional (projected) mass distribu- 
tion has units of mass (as one would expect). Therefore, the 
convergence and shear given by equations l|18p . I|19p and II20II 
are dimensionless. This is not surprising, and indeed this is 
true when one uses angular coordinates. The fiexion terms 
presented here are in dimensions of (angle) which is also 
consistent with those expressed in angular coordinates, how- 
ever the angle is in units of radians. This can be seen as the 
expressions for the first- and second-fiexion, equations l|37ll 
and l|4T|l respectively, are proportional to the ratio of the an- 
gular diameter distance from the observer to the lens plane, 
Dh, and the impact parameter, |^|. As discussed in section 
12.11 this ratio is . This implies that to compare the results 
presented in this pape r with those derive d using angular co- 
ordinates (for example lBacon et al.ll2006l l. one must convert 
from units of radians to arcseconds. As an alternative, in 
appendix [A] we show the equations for convergence, shear 
and fiexion expressed in angular coordinates. 



3 EXACT SOLUTIONS 

In this section we present analytic solutions of the equations 
expressed hitherto for various circularly symmetric matter 
distributions. For completeness, we first discuss the simplest 
case of a point mass, then move on to SIS, NFW and finally 
Sersic-law profiles. These exact solutions are then used in 
section |4] to investigate differences between the fiexion of 
various lenses, and we discuss how these fiexion terms can 
be used to constrain the density distributions. 



3.1 Point Mass 

Consider a Schwarzschild lens (i.e. a point mass) situated at 
the origin of the coordinate system on the lens plane. This 
implies that E = <5(|^|), where 8 is the Dirac delta function. 
Moreover, the projected mass distribution M = Ma is a 
constant for all |5| 7^ 0. The Einstein radius for such a system 
in angular coordinates is 



Be 



\ c2 



DsDl 



1/2 



(43) 



and the critical surface density and Einstein radius are re- 
lated according to 



Ms 



Ms_ 



(44) 



where we have defined £,e ■= DlOe to be the Einstein radius 
measured as a distance on the lens plane. 

The convergence is linearly proportional to the surface 
density of the system, equation HISII . implying it is triv- 
ially zero everywhere except |^| = 0. The shear, however, 
includes terms involving the projected mass of the system. 
Prom equations l|19p and l|20p the first and second compo- 
nents of the shear are 



71 



and 



72 



and the total shear is 



612 



2 

E 2i<h 

— e . 



(45) 



(46) 



Note that the final equation in l|46p is the usual expression 
for the shear associated with a point mass, however with 
an additional negative sign. This is a direct result of the 
calculation, however we note that the negative sign is due 
to the choice of coordinates. That is, rotating our coordinate 
system by ninety degrees implies the negative sign vanishes 
[i.e. (j) — > (j}+n/2 implies exp(2i(^) — exp{2i(j))]. Therefore, 
one is free to scale away the negative sign in equation ifiG]! 
by rotating the coordinate system, which yields the familiar 
result for the shear induced by a point mass lens. 

According to equations lf35|l and lf36|l. the components 
of the first-fiexion are the directional derivatives of the sur- 
face density, implying these vanish (for all |^| 7^ 0); 



J^i ^ ^ :f = 0. 



(47) 



As mentioned above, first-fiexion is associated with a shift of 
the centroid of the image with respect to the source. There- 
fore, = for a point mass implies the centroid of the 
image is unchanged. Second-fiexion, however, is associated 
with the "arciness" of the image, and one finds for the point 
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mass lens that 



Gi = 



—^2 (3^1 - 



Combining these, the total second-flexion is given by 



(48) 
(49) 

(50) 



Taking the magnitude of the above expression implies 
second-flexion decreases proportionally to the projected ra- 
dius cubed, \Q\ oc compared to the shear which 
decreases proportionally to the projected radius squared, 
ICr^. Therefore, as one gets further from the source, 
the shear term will come to dominate over the second- 
flexion, and this effect will become more negligible as one 
looks further from the origin (see flgure[2]). 

It is interesting to note that while the point mass has 
vanishing first-flexion and non-zero second-fiexion, it is rel- 
atively straightforward to also construct lens models where 
the second-flexion vanishes but the flrst-flexion is non-zero. 
Indeed by setting equation l|4T1l equal to zero, one can show 
that the second-flexion vanishes for models with E oc |^|^. 
Whilst this is obviously unphysical, as the surface density 
distribution increases monotonically as a function of radius, 
one still flnds that the flrst-flexion is non-zero everywhere 
(for all \i\ ^ 0) with \T\ oc j^j. 



3.2 Singular Isothermal Sphere 

Several observational studies of early-type galaxies using 
gravitational lensing have suggested that the total matter 
distribution (i.e. baryonic plus dark matt er) is well described 
by a nearly isoth ermal d ensity profile (Treu & Koopmana 



2002. 2004: Rusin et al 



. -- ^ . .20 03: Rusin & Kochanek 2 003 : 
Koopmaris "et al.ll2006l : ICzoske et al.l l2008l : IPve et all 12008 : 
Tu et al.ll2009ll. In deed t he strongest evidence for this is 



from lGavazzi et al.l |2007t l who reported weak lensing stud- 
ies of 22 early-type galaxies based on HST imaging, con- 
cluding that an isothermal distribution is consistent out 
to 100 effective radii. On larger scales, isothermal spheres 
have been found to fit the density profile o f galaxy clus- 
ters ijAthreva et all l2002l : lEttori et all l2002l : iKatgert etHI 
l2004h . although recent debate has suggested the NFW pro- 
file may provide a better fit (see section [373|l . Analytically, 
gravitational len sing studies have compar ed SIS and NFW 
profiles at first- jAVright fc Brainerdll2000l ) and higher-order 
ijSacon et al.l l2006l l. implying the SIS profile provides us 
with the simplest realistic profile to reproduce known re- 
sults. 

The SIS density profile is described in terms of the 
three-dimensional (deprojected) radius, r, as 



p(r) 



27rGr2 ' 



(51) 



where (J is the one-dimensional velocity dispersion. Project- 
ing equation II51II along the line-of-s ight gives the surface 
density (for example see lSinnev fc Tremaine 1987i ) 



(52) 



We note that the SIS profile, like the point mass, is singular 
at the origin. This implies that the following analysis is valid 
for all 1^1 7^ 0, however in reality one must quantify at what 
radius the system is a "weak lens" (see section|4|). By defining 
the Einstein radius in units of distance as S^e '■= DlOe, 
where is the standard Einstein radius for an SIS lens, one 
can show that the critical surface density, Ecr, is related to 
the Einstein radius as 



(53) 



By substituting the above equations into ifTSll . we can show 
that the convergence falls-off linearly with the distance from 
the center of the coordinate system; 

Moreover, the components of the shear have a similar be- 
haviour with the extra angular dependence 



71 



72 = 



(55) 



(56) 



As we know the shear is a spin-two field, we are essentially 
only interested in the strength of the respective fields. There- 
fore, for the remainder of the article, unless explicitly stated, 
we shall just be working with the magnitudes of these fields. 
Now, as 7 = I7I exp{2i(f>), one can show from the above that 

I £,E _ ds 

^ 2l^ ^ 26' 
The first-flexion components are 

-Dl£.e . 



(57) 



^ -D^E 



and 



^2 



'<2. 



(58) 



First-flexion is a spin-one field, implying = exp(j0). 
Equations l(58|l imply the first-fiexion falls-off proportionally 
to the distance squared; 

Dl£,e ^Oe 



2|5|2 2612 • 
The second-fiexion components are 
SDl^e . 



Q2 = 



?.Dl^e 



6 (G - 36') , 
6 (3C? - &) , 



(59) 



(60) 



(61) 



2Giei 



which is a spin-three field, G = \G\ exp{Si(j}) , implying 
second-flexion also decreases proportionally to the radius 
squared 

1^1 - = (62) 

All of the above equations that are written in angular coor- 
dinate s are consistent with those presented in iBacon et al.l 
Hooi). These analytic expressions are plotted for a specific 
mass in figure [2] where they are compared with the NFW 
and Sersic-law profiles. 

3.3 Navarro-Frenk- White Profile 

Compared to the SIS profile, the NFW profile introduces an 
extra parameter into the scaling of density distributions. The 
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concentration, c, which is defined as the ratio of the three- 
dimensional virial radius to the three-dimensional scale ra- 
dius, c = VA/rs, is a function of the particular cosmology 
being used. In section [4] we compare the NFW profile to 
the other profiles being analysed here, as well as compar- 
ing the effects of varying the concentration on the first- and 
higher-order lensing phenomena. The NFW profile is gener- 
ally given in terms of the three-dimensional radii, r. 



p{r) 



5c Pc 



(r/rs) (1 + r/rsf ' 
where pc is the critical density of the Universe and 

5 

3 ln(l + c)-c/(l + c)' 



(63) 



(64) 



Projecting this ont o the two-dimensi onal radius, |^|, gives 
the surface density ljBartelmannlll99(j l 



E = 



2pc5crl 



[1-H(|?|)] 



where we have defined the following function 
= arctanh4/ 



H(ICI) 



:arctan 



lei- 



(65) 



(66) 



Integrating the surface density gives the projected mass dis- 
tribution 



M = A-JvpcScr's 



ln|L+H(|,|: 



(67) 



The convergence and shear can now be expressed simply 
as functions of the above expressions. From equation l|2ip . 
the convergence is simply k = E/Scr, where the surface den- 
sity is given by equation l|65p . The total shear is determined 
by equation l|22p as 



l7l 



2pc5crl 



Sc.(iei= 



H-2 1 



2rs 



,(68) 



where H = H(j^j). Equation l|37p implies first-fiexion is found 
by differentiating the surface density, which can be shown 
to be 



1^1 = 



^2DLPc5cri 



T.cr\i\m'-ri 



(69) 



and equation l|4ip implies the second-fiexion for the NFW 
profile is given by the expression 



\Q\ = 



2DLPc5cri 



E„.|C|(|C|2-ri) 



2^2 



8 1 



+31 



■2\i\ 



2Qrl + 8 



(70) 



Despite the mass of the system being infinite, one can show 
that the convergence, shear, first- and second-fiexion all tend 
to zero as |^| — > oo. These profiles are plotted against the 
SIS and Sersic-law profiles in figure El and we also look at 
the dependence of the concentration in figure [3] and section 



3.4 Sersic Profile 



It has long been argued that a Sersic-law ijSersid Il968l ) 
provides a remarkably good fit to luminosity profiles of 
early-type galaxies, ranging in size f rom dwarf galax- 
ies to the lar g est elliptical galaxi e s JCaon et al. _ 19931: 
GrahamI [iooH: iGraham fc Guzma3 l2003l : I Graham et al.l 



20031 : iTruiillo et al.ll2004l'l. For a concise reference to Sersic 



quantities, see I Graham fc Driverl l|2005l 'l. Recently, a Sersic- 
law has also been shown to provide a good fit to three- 
dimensional density profiles jNavarro et al.ll2004|). and a lso 
to projected surface density profiles I Merritt et al.l ioOSll of 
dark matter halos. In a series of papers jlvlerritt et al.l 20061 : 
iGraham et"aDl2006al lbh it has further been shown that pro- 
jected Sersic surface density profiles provide the best fit to 
sim ulated galaxy- a nd cluster-sized dark matter halos. 

ICardond l| 20041 ) first analysed the Sersic profile in the 
gravitational lensing context, showing that mass estimates 
using lens reconstruc tions is highly dependent on the choice 
of Sersic parameter. lEHasdottir fc Mollerl HoQT') compared 
gravitational lensing for Sersic and NFW profiles, and found 
that mass estimates may differ by up to a factor of two, de- 
pendent on the choice of density profile and Sersic index. 
In the weak lensing regime, they did this by looking at the 
shear of both profiles. We take this a step further by also 
analysing higher-order lensing terms. In this section, we pro- 
vide the first explicit representation of fiexion terms for the 
Sersic-law profile, and in section|4]we compare these results 
to those of the NFW and SIS profiles. 

The Sersic profile is defined in terms of the surface den- 
sity 



In 



1/n 



(71) 



where Ee is the surface density at the effective radius, ^e. 
The constant n is the Sersic shape parameter which de- 
scribes the shape of the profile and 6„ is a function of n 
that is chosen such that the effective radius contains half of 
the projected mass of the system. Analytically, this is given 
as the solution of r(2n) — 2V{2n, bn), where 



r(Q, x) 



(72) 



is the lower incomplete gamma function and r(a) := 
lim^j^oo r(a, x) is the complete gamma function. As such, 
bn can be reasonably approximate d to bn = 2n — 1/3 - I- 
4/(405n)-HC»(n"2) for 0.5 < n< 10 |Ciotti fc Bertinlll999l l. 
Integrating the surface density gives the projected mass 



M = 27rn^E4,"r(2n, Z) 



(73) 



where Z = 6„ (|C!/Ce)'/". One can see that the Sersic-law 
has one more parameter than the NFW profile, a point we 
discuss in more detail in section [4] and also appendix [BI 

After much algebra, one can show that the magnitudes 
of the convergence and shear for the Sersic profile can be 
expressed as 



■ exp < bn 



l/n 



(74) 
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h\ = ^ exp <( h 




(75) 



Furthermore, the magnitudes of the first- and second-fiexion 
terms are 



1^1 



1/n 



exp < bn 



l/n 



(76) 



6n neh'/" 



X exp < bn 



l/r. 



(77) 



4 PROFILE COMPARISONS 

We compare gravitational lensing effects of the various 
density distributions by holding the virial mass of each 
system constant, which is the same method used by 
IWright fc BrainerdI l|200(]h for comparing first-order lens- 
ing properties of the NFW and SIS profiles. The three- 
dimensional (i.e. deprojected) virial radius, ta, is defined as 
the radius inside which the average density of the halo is A 
times the critical density of the Universe, implying the three- 
dimensional virial mass is Ma = ^TvApcr^/S (throughout 
the remainder of the article we use A = 200). Constructing 
the lens models then requires the three-dimensional (depro- 
jected) density distribution for each profile, as well as the 
three-dimensional mass distribution. In general, these two 
equations can be inverted to find va, and also the various 
parameters associated with the individual profiles (for ex- 
ample the velocity dispersion, a, for the SIS profile). The 
Sersic profile is a little more difficult to treat with this pro- 
cedure as analytic forms of the deprojected density and mass 
functions do not e xist. As such, we use the a nalytic approx- 
imations given bylPrugniel & Simien' fl99^, as well as em- 
pirical relations from [Graham et al. (2006l3). Details of the 
way in which we construct Sersic-law density distributions 
are provided i n appendix [bI 

Following iBacon et anij2006l l. we use a fiat ACDM cos- 
mology with Clrn = 0.3, f^A = 0.7 and h = 0.72. We place 
the lensing mass at a redshift of zl = 0.35 and the source 
at zs = 0.8 as these values correspond to Dls/Ds — 0.5. 
To compare the four density profiles, we use a lens mass of 
M200 = lO^^/i"^M0. For the NFW profile, the concentra- 
tion factor, c, which is the ratio of the virial radius to the 
scale radius (see section [33|l . is a function of the cosmology 
and the redshift of the lens which, for the above system, is 
evaluated to be c = 7.20. For the Sersic model we use empir- 
ically derived relations between the Sersic shape parameter , 
n, and the mass of the system: from I Graham et sl\ lj2006lJ ). 
their equation (12), we find the Sersic shape parameter for 
a galaxy of mass M200 = 10^^ Mq is n ~ 8.6. 

Figure [2] shows the lensing properties of the 
Schwarzschild lens (dotted red line), SIS (dashed green line), 
NFW (blue dashed-dot line) and Sersic (thick black line) 



mass distributions. The Schwarzschild lens has an Einstein 
radius of — 11.5 kpc, which implies that for impact pa- 
rameters inside this radius the results presented are not in 
the weak lensing regime. The most striking feature in these 
plots is that at large distances, the convergence and first- 
fiexion for the SIS are significantly larger than for the NFW 
and Sersic profiles. Bearing in mind that the convergence 
is linearly proportional to the surface density, this implies 
that at large distances the surface density of the SIS is also 
significantly larger than the NFW and Sersic profile. This is 
consistent with the fact that the outer-logarithmic slope for 
NFW profiles is —3, while the slope of an isothermal sphere 
is —2. Additionally, the slope of a Sersic-law profile depends 
on the specific shape parameter, n, and varies as a function 
of the distance from the centre of the profile, i.e. Inter- 
estingly, whilst these features are evident in the convergence 
and the first-fiexion, they are less apparent in the shear and 
second-fiexion. 

It is worth exploring the extent to which the above 
properties of the NFW and Sersic profiles depend on the 
concentration and shape parameters respectively. The con- 
centration parameter in the NFW profile, c, is defined as 
the ratio of the virial radius to the scale radius, which is a 
function of the specific cosmology. Figure [3] shows the effect 
of a varying concentration on the different lensing proper- 
ties. A lensing mass of M200 = 10^°h~^ Mq is used, with 
values of c = 4, 8, 12, 16, 20, 24. For a Schwarzschild lens of 
this mass at these distances the Einstein radius is — 1.2 
kpc, which again gives us a scale on which a weak lensing 
treatment is appropriate. 

It is apparent from figure [3] that the lensing properties 
of an NFW profile are not linearly effected by the concen- 
tration parameter. Indeed as c becomes larger, the lensing 
properties begin to converge, implying that at large c the 
plots become indistinguishable. Moreover, the effect of the 
concentration parameter is greater at small distances from 
the centre of mass of the lensing galaxy. This is seen most 
pertinently in the convergence where at j^j ~ W^e the lines 
are indistinct. It is also interesting to note that variations 
in the concentration parameter cause the shear and second- 
flexion to change signiflcantly more than the convergence 
and first-fiexion. Therefore, given a series of lensed images, 
and assuming an NFW fit to the density profile, one can 
learn more about the specifics of the profile from the shear 
and the second-fiexion than from the convergence and the 
first-fiexion. There is likely to be a degeneracy between the 
mass of the lensing galaxy and the concentration parameter 
if there are only a handful of images. Our results imply that 
using both the shear and the second-fiexion may be able to 
break this degeneracy so c and M can be obtained. 

Figure [4] shows the lensing properties of Sersic profiles 
for various values of the Sersic shape parameter, n. The lens- 
ing mass is again M200 ~ lO^^/i^^M©, implying the point 
mass lens has ^e — 1.2 kpc, and we vary 1 < n ^ 9. The 
striking feature of these plots is that the shear and second- 
fiexion have significantly less dependence on the shape, n, of 
the Sersic profile than the convergence and the first-fiexion. 
Indeed, zooming in on the vertical scale by an order of mag- 
nitude for the shear and second-fiexion plots reveals that the 
dependence on n can only be seen for |^| ^ ^e- Moreover, 
the effect of the shape parameter on the convergence and 
first-fiexion increases as one moves further from the source, 
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Figure 2. Convergence, shear, first- and second-flexion for Sciiwarzsctiild lens (dotted red line), SIS (dashed green line), NFW (dash-dot 
blue line) and Sersic profile (thick black line). The three-dimensional virial masses of each of these profiles is 10^^ Mq, with the lens 
and source placed at zl = 0.35 and zs = 0.8 respectively. As an indicator of the weak lensing regime in these plots, the Einstein radius 
for the Schwarzschild lens is — 11-5 kpc 



although this may be difficult to detect observationally due 
to the relative size of the signal being significantly weaker at 
large separations from the lensing mass. The fact that there 
is no dependence on the specific shape of the Sersic profile 
on the shear and second-flexion implies that these two prop- 
erties can be used to derive the mass of the lensing object, 
whilst the convergence and first-fiexion can then be used to 
derive n. Weak gravitational lensing thus provides an inde- 
pendent method for deriving the masses of Sersic-law galax- 
ies and clusters which only weakly depends on the specific 
shape of the profile, provided the correct lensing properties 
are utilised, i.e. the convergence and second-fiexion. 



5 CONCLUSION 

We have derived general equations governing the conver- 
gence, shear, first- and second-fiexion for circularly sym- 
metric gravitational lenses in terms of the surface density 
and projected mass of the lens. We have shown that the 
components of the first-fiexion are simply the directional 
derivatives of the surface density, while the second-fiexion is 
a slightly more complicated function of the surface density. 



its gradient and also the projected mass distribution. By 
applying the formalism to specific lens models, in particu- 
lar a Schwarzschild lens, a singular isothermal sphere (SIS), 
Navarro-Frenk- White (NFW) profile and a Sersic profile, we 
have compared the signature each profile has on each of the 
lensing terms as a function of the radial impact parameter. 
Whilst the NFW and SIS pr ofiles have been conipared pre- 
viously, both a t linear-order dWright fc Brainerdlfioool ) and 
also for fiexion ijBacon et al.ll2006l l. to the best of our knowl- 
edge this is the first time fiexion for Sersic-law profiles have 
been presented. In particular, we showed that the shear and 
second-fiexion effects for a Sersic profile are systematically 
larger than for the NFW and SIS profiles (figure El . This 
implies that one must be careful about the specific mat- 
ter distribution assumed when attempting to analytically 
reconstruct the total mass and mass profile of a circularly 
symmetric lens. 

An SIS profile is uniquely determined by designating 
the mass at a certain radius (e.g. virial mass). However, 
both the NFW and Sersic models have extra parameters 
that determine their exact shape. The concentration of the 
NFW profile is the ratio of the virial radius to the scale ra- 
dius, which is a function of the particular cosmology being 
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Figure 3. Convergence, shear, first- and second-fiexion for NFW profiles with ranging values of the concentration, c, with M200 = 
W^'^h~^ Mq. Note that a Schwarzschild lens of this mass at this distance has an Einstein radius of £,e — 1-2 kpc, implying the weak 
lensing regime is somewhat beyond this limit. The values of the concentration shown are c = 4, 8, 12, 16, 20, 24. As the concentration 
increases, the curves converge on one another. The shear and the second-flexion show distinctions between the various values of the 
concentration better than the convergence and first-flexion. 



considered. We showed that the convergence and first-flexion 
are relatively weakly dependent on the concentration param- 
eter used, compared with the shear and the second-flexion 
(see flgurelS}. The Sersic proflle is parametrized by the Ser- 
sic shape parameter, n. Counter to the NFW concentration 
parameter, the Sersic shape parameter is more heavily de- 
pendent on the convergence and flrst-flexion, whereas the 
shear and second-flexion vary little as n changes (flgure|4|). 
These properties therefore provide the opportunity to di- 
rectly measure both the mass of lensing galaxies as well as 
the speciflc shape of their density proflles. Given a limited 
supply of images around a speciflc galaxy, higher-order lens- 
ing terms may provide the ability to give extra constraints 
on the individual proflles of the galaxies. 

The gravitational lensing systems explored in this ar- 
ticle are idealized in that the projection of their mass dis- 
tributions are circularly symmetric. Moving beyond circular 
symmetry requires the numerical solution of the thin-lens 
gravitational lens equation. One is then free to study the ef- 
fect of non-circular lens models by including anisotropics in 
the lensing galaxy, and also to look at the effect the size of 
the source has on the shape of the flnal image. Whilst these 



applications will be extremely useful for determining the 
mass distributions of various galaxies using next-generation 
gravitational lensing surveys, the higher-order gravitational 
lensing effects will also be useful in studying the dark matter 
and dark energy content of the Universe. Moreover, although 
the study of flrst-order gravitational lensing through N-body 
simulations has been aro und for some time (see for exam- 
ple the recent review bv lMunshi et al.]|2008l l. the study of 
flexion in these systems has scarcely been broached. In this 
way, one can study the expected probability distributions 
for flrst- and second-flexion as a function of the speciflc cos- 
mology. The advantage of using flexion as opposed to linear 
gravitational lensing effects is that a systematic bias is taken 
out of the study as only one reasonable assumption about 
the source object is required - that it is not intrinsically 
flexed. 
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Figure 4. Convergence, shear, first- and second-flexion for Sersic profiles with ranging values of the Sersic shape parameter, n, with 
M200 = IO^O/i-IMq. As with fi gure[3l a Schwarzschild lens of this mass at this distance would have an Einstein radius of — 1-2 kpc, 
implying the weak lensing regime is somewhat beyond this limit. The shear and second-fiexion vary extremely weakly with a change in 
n, whilst the convergence and first-flexion are more heavily dependent on this parameter. 



APPENDIX A: CONVERGENCE, SHEAR AND 
FLEXION IN ANGULAR COORDINATES 

It is instructive to show the main equations from sections 
1 2. 2 1 and 1 2. 3 1 expressed in angular coordinates, {0i, 62), where 
= DhOi. From the definition of the projected mass func- 
tion in terms of the surface density, equation implies 
that the projected mass expressed in angular coordinates 
is related to the mass in distance coordinates according to 
M(|5|) = DlM{\e\), where \e\ = ^Of+Ol It is trivial to 
show that the expression for the convergence, equation l|18p 
goes unchanged, however the components of the shear, equa- 
tions (fT9ll and ll20ll . written in angular coordinates become 



71(^1,^2) = 
72(6*1, 6)2) = 



7rE„.|e|4 
26)i6>2 



M) 
M) 



(Al) 
(A2) 



where for the remainder of this appendix E = E(|0|) and 
M = M{\e\). The magnitude of the total shear is 



17(^1,^2)1 



1 



UEier 



M) 



(A3) 



The components of the first-flexion, equations l|35p and l|36p. 
are given in angular coordinates as 

1 9E 
Ed- d9^ ' 
1 dT. 

Tier 002 ' 

implying the magnitude of the first-flexion is 
1 dT. 



^1(91,62) = 
^^2(^1,^2) = 



(A4) 
(A5) 



2)1 = 



Ecr d\e\ 



(A6) 



Where previously the flexion terms were expressed in units 
of distance, one can now see first-flexion has units of 
(angle) It is trivial to see that the following three ex- 
pressions governing the second-flexion in terms of angular 
coordinates are also expressed in these same units. 



Gi {61,62) 
Q2 {61,62) 
\Q{6i,62)\ 



di {6l ~ 361) 
7rE,,|e|e 

62 {36f - el) 



7rEcr|6l|e 
1 / dE 

E„. V dl^l ' 



_ dT, 

'd\e\ 

_ dT 

diei 



4E 4A/ 



47rEi6'r +4j|A7,) 
47rE|6'l^ +44A8,) 
(A9) 
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APPENDIX B: CREATING GALAXIES WITH 
SERSIC PROFILES 

The Sersic profile is defined in terms of the projected sur- 
face density, equation l|7ip . However, to compare the lens- 
ing effects of the Sersic profile with other density distri- 
butions, one requires the deprojected form of the density 
and mass distributions (see the discussion in section |4]). 
Whilst an analytic form of the deprojected Sersic profile is 
not available, an analytic approximation has been provided 
b;jPrugniel fc Simien (199 7), an d has further been explo red 
in d etail bvlMerritt et all ll2006l l: ICraham et aP lj2006al lbh. 

IPrugniel &: Simie nl l|l997l ) showed that the three- 
dimensional density distribution associated with the surface 
density of the Sersic profile, given by equation l|7ip . can be 
approximated as 



p(r) = 



exp < —bn 



1/n 



- 1 



(Bl) 



Here, pe is the three-dimensional density at the effective 
(projected) radius ^e- The function p = p{n) is utilised 
to ensure that the projection of equation l|Bip relates as 
closely as possible to the projected Sersic profile, i.e. equa- 
tion Il7ip. for the range 0.6 ^ n ^ 10. This was first given by 
iLima Neto et al.l ((1999) as p = 1. - 0.6097/n -f 0.054 63/w^ , 
and a goodness of fit is shown in lMerritt et al.l l|2006l ). Inte- 
grating the three-dimensional density distribution over the 
volume gives the three-dimensional mass distribution, 



MiD{r) 



A ^3 b„ 



-(3-p)nT 



(3 - p) n, 6„ ( — 



l/r. 



.(B2) 



ICardonj l|2004l ) first looked at gravitational lensing for 
a Sersic profile. He discussed the need to reduce the pa- 
rameter space of the system in order to build a graviational 
lens from a Sersic model. To that end, he used an empiri- 
cal relation governing the deprojected effective radius, the 
central surface brightness and t he Sersic shape parameter, 
n (see equation A. 4 and A. 5 of [Cardofll (20041 )). It is pos- 
sible for us to also use this relation, and subsequently con- 
vert the central surface brightness into a density by invoking 
more empirical relations and also assuming a mass-to-light 
ratio. Essentially, this procedure has already been completed 
for the Prugniel & Simien model bv lGraham et al.l l|2nnfibh 
[their equations (13) and (14)]; 



logic Pe^k- 2.5 login 



(B3) 



Here, is in units of kiloparsecs, pe is in solar masses per 
cubic parsec and fc is a constant which is 0.5 for luminous 
elliptical galaxies and galaxy-sized dark matter halos (with 
log Re > 0.5) and 2.5 for cluster-sized dark matter halos 
(with logi?e > 1.5). Finally, the three-dimensional density 
at the effective (projected) radius is related to the two- 
dimensional surface density at by 

r(2n) 



2C.r[(3-p)n]- 

As mentioned in section [H to compare profiles we spec- 
ify the virial mass of the system. Ma, implying we know 
the virial radius, ta. Substituting this into equation ljB2p . 
together with equation l|B3P implies we have an equation 
for 5e as a function of the Sersic shape parameter, n. This 
equation is not analytically invertible due to the presence 



of the incomplete gamma function, however it can be solved 
numerically for given values of n. Therefore, once this equa- 
tion is solved, we know and hence Pe, which can both 
be substituted into equation l|B4p to give Se. The projected 
Sersic profile l|7ip can finally be evaluated, along with its 
various derivatives and also the two-dimensional mass, im- 
plying all of the lensing quantities can be evaluated. 
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